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INTRODUCTION 


The  objective,  scope,  and  technical  effort  for  this  program,  as  stated  in  the  contract 
STATEMENT  OF  WORK,  were: 

1.0  Objective.  The  objective  of  this  work  is  to  develop  a  theoretical  description  for 
nonlinear  guided-wave  degenerate  four-wave  mixing. 

2.0  Scope.  The  effort  will  consist  of  theoretical  investigations  of  the  nonlinear 
properties  of  nonlinear  waveguides.  The  studies  will  include  the  investigation  of 
four-wave  mixing. 

4.0  Technical  Effort.  The  contractor  shall  review  the  existing  theory  for  nonlinear 
waveguides  including  solilon-based  nonlinear  guided-wave  phenomena.  The  contractor 
shall  develop  a  theoretical  description  and  assessment  for  degenerate  four-wave 
mixing  in  nonlinear  waveguides. 

The  first  phase  of  the  program  consisted  of  reviewing  the  literature  on  the  two 
problems  of  interest  to  this  contract,  namely 

1.  degenerate  four-wave  mixing  in  waveguides. 

2.  soliton  transfer  between  waveguides. 

The  second  phase  of  the  program  consisted  of  original  theoretical  calculations  on 
both  of  the  above  problems. 

REVIEW  OF  DEGENERATE  FOUR-WAVE  MIXING  IN  WAVEGUIDES 

Degenerate  four-wave  mixing  (DFWM)  in  fiber  and  planar  waveguides  has  been 
discussed  theoretically  and  demonstrated  experimentally.  Previous  work  in  this  area, 
up  to  early  1986.  is  summarized  in  Ref.  1.  In  the  case  of  fibers,  all  of  the  truly 
original  work  was  done  before  1983.  This  waveguide  system  is  of  marginal  interest 
since  the  interacting  beams  must  travel  along  the  fiber  axis;  that  is.  fibers  correspond 
to  a  one-dimensional  system.  The  possibilities  for  working  in  planar  waveguides  are 
much  more  diverse  because  the  interacting  beams  can  either  lie  in  the  plane  of  the 
film,  or  be  incident  onto  it  from  above  or  below.  Therefore,  in  this  repr.t,  we 
concentrate  on  degenerate  four-wave  mixing  in  planar  waveguide  structures. 

There  are  essentially  three  viable  geometries  for  performing  degenerate  four- 
wave  mixing  in  planar  waveguides: 

1.  All  four  waves  guided. 

2.  Two  externally  incident  pump  beams,  and  guided  probe  and  conjugate 
beams. 
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3.  Two  guided  pump  beams,  and  external  probe  and  signal  beams. 

The  first  case,  all  four  waves  guided,  has  been  analyzed  by  at  least  two  groups, 
ourselves  included.1'*  In  fact,  two  experiments  have  been  reported.4'*  The  first  used 
the  nonlinear  material  carbon  disulphide  as  the  cover  material  above  the  guiding  film. 
Because  the  fraction  of  power  guided  in  the  cover  medium  is  small,  the  efficiency 
for  this  geometry  was  low.  In  the  second  experiment,  the  waveguide  itself  consisted 
of  the  nonlinear  medium,  namely  semiconductor-doped  glass,  and  the  efficiencies  were 
very  good.  For  peak  powers  of  tens  of  watts,  efficiencies  of  10%  have  been 
achieved  with  picosecond  pulses.  It  is  because  the  guided  beams  are  strongly 
confined  in  one  dimension  that  the  reflectivities  per  unit  input  power  for  the  pump 
beams  are  large.  However,  for  a  guided  wave,  the  phase  information  can  be  coded 
in  only  one  dimension;  that  is,  along  the  wavefront,  perpendicular  to  the  propagation 
direction. 

The  geometry  with  two  external  pump  waves,  incident  from  above  and  below 
the  guiding  film,  and  a  guided  probe  beam  leading  to  a  guided -wave  signal  beam  has 
also  been  treated  theoretically.6  However,  the  applications  of  such  a  geometry  are  not 
obvious. 

The  geometry  which  best  utilizes  guided  waves  for  phase  conjugation  is  the 
third.  As  illustrated  in  Fig.  1.  the  contra-propagating  pump  beams  are  guided, 
meaning  that  small  total  powers  are  required  to  obtain  high  waveguide  intensities. 
The  probe  beam  is  incident  from  above  the  waveguide  and  a  phase  conjugate  signal 
is  obtained  both  on  reflection  and  transmission  through  the  waveguide.  This 
particular  case  has  been  neither  analyzed  nor  demonstrated  experimentally.  This  case 
is  treated  theoretically  in  the  next  section. 

DEGENERATE  FOUR-WAVE  MIXING  WITH  GUIDED  PUMP  BEAMS 

The  geometry  shown  in  Fig.  1  was  analyzed  in  the  weak  conjugate  signal  limit. 
The  counter-propagating  pump  beams  are  both  assumed  to  be  TEq  guided  waves;  that 
is.  the  electric  field  polarization  is  the  y-direction.  (We  do  not  treat  the  TM  case 
here,  although  we  do  not  expect  anything  to  be  substantially  different  other  than  the 
existence  of  a  Brewster  angle  phenomenon.)  The  probe  signal  is  incident  from  above 
the  thin-film  waveguide  in  the  form  of  a  plane  wave  at  an  angle  8e  relative  to  the 
normal.  In  addition  to  the  usual  reflected  DFWM  signal,  there  is  also  a  "transmitted" 
DFWM  signal  attributable  to  multi-reflections  inside  the  film. 
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Fig.  1.  The  DFWM  geometry  analyzed  here.  The  pump  beams  are  guided  by  a  film  of 
thickness  h.  A  plane-wave  beam  is  incident  from  the  cladding  side  and  DFWM  is  obtained 
both  on  reflection  and  transmission  through  the  film. 


The  calculation  of  the  nonlinear  polarization  associated  with  the  DFWM  signal  is 
relatively  straightforward.7  Including  the  degeneracy  factor. 

PyNL(r)  .  4<0  02E(x)  E+ME_ME,nc‘M  .  (1) 

where  n^  is  defined  in  the  usual  way  as  a  field-dependent  refractive  index  of  the 
form  n(|E|)  -  n„  +  n^l  E| 1  (and  E  is  the  local  field).  Here  E„M  and  Ejw)  are  the 
counter-propagating  pump  beams,  and  E,nc(u))  is  the  plane-wave  incident  field. 

The  strength  of  the  DFWM  signal  on  reflection  and  transmission  is  calculated  by 
solving  the  polarization-driven  wave  equation  for  fields  with  wavevector  component 
-Kpk„  parallel  to  the  surface;7  i.e.. 


1  d* 
dx1 


Ey(W.X) 
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<o  y 


where 

E.1,.0  -  \  £,<«.*)  «i(-,Pk'2  -  "*l  *  c.c  . 

The  form  of  the  DFWM  signal  fields  far  from  the  film  is  given  by 

nc:  Ey(w.x)  -  Dce”‘*c,CoX  when  x  <  0  . 
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3 


(2b) 


a  Fop 

tkl 


tton_ 
(3a)  .  - 


(3b)  U°n/ - 

lllty  Cod»» 
.nitTand/or 
Plat  I  Speeial 

l.  i  1 


□  □ 


Noting  that  the  intensity  S  is  given  by  0.5no«o|  Ey| 1  and  the  (intensity)  DFWM 
reflection  and  transmission  coefficients  are  given  by 

R-  it--  |Dc|*-nrP,P.  .  (4a) 

°it\c 

Tm  r~m  •  (4b) 

“inc  uc 

The  details  of  the  solutions  can  be  found  in  the  appended  paper.7 

Numerical  calculations  were  performed  on  two  representative  guided-wave  film 
geometries,  both  utilizing  typical  nonlinear  organic  parameters.  In  the  first  case,  the 
following  parameters  were  assumed:  nc  «  1.0.  nr  -  1.6.  ng  -  1.45,  X  -  1.06  jim. 
n2ir  ■  10~16  to2/ W  and  n2Jt  -  0.0  In  the  second  case,  the  effect  of  a  nonlinearity  in 
the  substrate  instead  of  the  film  was  examined,  with  n2if  -  0.0  and 
n2h  -  10-“  m*/W.  Relatively  large  oscillations  in  the  DFWM  cross  section  were 
observed  for  both  cases,  and  two  additional  calculations  were  performed  with  a 
smaller  difference  between  the  film  and  substrate  refractive  indices,  namely  nr  -  1.6 
and  n,  -  1.59. 

Figures  2  and  3  show  ijr  and  i>t,  respectively,  as  a  function  of  incidence  angle 
6C.  For  film  thicknesses  near  waveguide  cutoff,  the  externally  incident  beam  does 
not  show  any  resonances  across  the  film  and  the  variation  in  the  cross  sections  T)r 
and  i \  is  smooth  with  angle.  However,  as  the  film  thickness  is  increased  and 
resonances  begin  to  appear  across  the  film,  distinct  maxima  and  minima  in  the  cross 
section  appear,  their  number  increasing  with  film  thickness.  Since  the  transmitted 
DFWM  signal  requires  reflection  at  the  lower  film  boundary,  the  oscillations  are 
much  larger  on  transmission  than  on  reflection.  Therefore  it  appears  feasible  to 
optimize  the  signal  for  a  given  film  thickness  by  varying  the  incidence  angle. 

The  variation  in  t)  with  film  thickness  at  a  45°  incident  angle  is  shown  in 
Fig.  4.  As  expected,  oscillations  again  occur  because  of  multi-reflections  of  both  the 
incident  and  signal  beams  in  the  film.  When  the  index  difference  between  the  film 
and  substrate  is  reduced,  the  amplitude  of  the  oscillations  is  also  reduced,  as  shown 
in  Fig.  5. 

The  effectiveness  of  using  a  linear  film  on  a  nonlinear  substrate  is  examined  in 
Fig.  6.  Clearly  the  cross  section  decreases  rapidly  from  its  cutoff  value  and  falls 
orders  of  magnitude  below  that  obtained  on  reflection. 
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Fig.  2.  The  reflection  DFWM  cross  section  tir  versus  angle  of  inidence  0, 
different  film  thicknesses.  Here  nc  *  1.0,  nj  *  1.6,  ns  »  1.45,  n2Jf  *  1x10' 

n21s  *  0  0 


Fig.  3.  The  transmission  DFWM  cross  section  r}t  versus  angle  of  inidence  0, 
different  film  thicknesses.  Here  nc  *  1.0,  n/  *  1.6,  ns  *  1.45,  n2/f  *  IxlO' 16 

n21s  =  00 


for  three 
n?!W  and 


for  three 
rrPIW  and 


5 


MMCTtON 


Fig.  4.  The  reflection  (r\r)  and  transmission  (i\t)  DFWM  cross  sections  versus  film 
thickness  for  6C  *  45°.  Here  nc  «  1.0.  nt  *  1.6.  ns  *  1.45,  njit  *  1x10 -,s  rtlW  and 
n21s  *  0.0  " 


Fig.  5.  The  reflection  ( i)r )  DFWM  cross  section  versus  film  thickness  for  6C  *  45°.  Here 
nc  =  1.0,  nj  *  1.6,  ns  «  1.59,  n2ff  *  lxlO‘H  n?/W  and  n2/s  *  00 


Fig.  6.  The  reflection  fr)r  s  R)  and  transmission  (t)t  b  T)  DFWM  cross  sections  versus 
film  thickness  for  &c  *  45°.  For  the  solid  line.  nc  *  1.0,  nj  *  1,6.  ns  ~  1.45, 
n2/ s  =  1x10' 16  nP/W  and  n/n  *  0.0.  For  the  dotted  line,  nc  *  1.0,  ni  =  1.6,  n5  *  1.59, 
n21s  “  M0-"  nP/W  and  n///  *  0.0 


We  now  compare  the  detailed  behavior  shown  in  the  preceding  figures  with  an 
approximation  based  on  plane-wave  analysis.  It  is  well  known  that  for  plane  waves 
in  the  small  signal  limit,* 


R  -  -  16t*^n|i  S+S.  ,  (5) 

where  n  -  n<,  +  n2iS  and  L  is  the  interaction  distance.  To  a  useful  approximation 
S.  -  P./hlff,  where  heff  is  the  effective  waveguide  thickness.  Therefore,  rewriting 
for  the  guided- wave  case. 


2 

R  -  16ir>  p+p  .  (6) 

n«rr  K 

should  provide  a  useful  approximation.  For  thick  films,  heff  #  h  *  L  and  the  cross 
section  coefficient  becomes  independent  of  film  thickness.  Estimating  for  the  present 
case,  we  obtained  0.16xl0'17  for  the  cross  section,  in  good  agreement  with  the 
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"average"  value  of  the  oscillations  of  0.2x1 0_IT  shown  in  Figs.  4  and  5.  Despite  the 
fact  that  h  ■  L  is  usually  good  to  at  best  ±  25%,  Eq.  (6)  does  give  a  useful  value  for 
the  cross  section.  In  addition,  it  predicts  the  asymptotic  behavior  with  increasing 
film  thickness.  Furthermore,  since  h#ff  >1  cutoff,  where  the  substrate  field 
degenerates  into  a  plane  wave  traveling  parallel  to  the  surface.  R  ■*  0  at  cutoff  can 
also  be  understood.  Finally,  it  is  clear  that  the  asymptotic  value  of  the  cross  section 
with  film  thickness  depends  primarily  on  the  nonlinearity  of  the  film,  and  not  the 
details  of  the  waveguiding  structure. 

Equations  5  and  6  also  provide  some  insight  into  the  maximum  values  of 
reflectivity  that  might  be  available.  Assuming  power  densities  approaching  damage 
values,  for  example  10  GW/cm!,  and  I -jim-thick  films,  this  corresponds  to 
P  a  107  W/m,  leading  to  a  maximum  reflectivity  of  0.02  (2%).  (For  a  1-cm  guided- 
wave  beam,  a  peak  power  of  only  100  KW  is  required  for  the  pump  beams.)  Hence, 
for  this  material  system,  pump- wave  depletion  can  be  ignored.  For  semiconductor 
materials  (which  also  exhibit  large  absorption),  reflectivities  in  excess  of  100%  should 
be  possible. 

COMMUNICATION  BETWEEN  WAVEGUIDES  VIA  SOLITON  TRANSFER 

One  way  to  communicate  information  and  enforce  phase  locking  between  two 
waveguides  is  via  the  transfer  of  energy  between  them.  The  usual  way  to 
accomplish  this  is  to  place  two  parallel  waveguides  in  close  proximity  so  that  their 
fields  overlap  and  energy  exchange  occurs  with  propagation  distance.  In  the 
integrated  optics  community,  such  an  arrangement  is  called  the  "directional  coupler." 
Another,  novel  approach  is  to  use  the  exchange  of  spatial  solitons  between  the  two 
waveguides.  It  is  this  approach  which  is  examined  here. 

Interaction  between  waveguides  via  the  exchange  of  spatial  solitons  only  became 
a  possibility  about  two  years  ago  when  it  was  shown  that  spatial  solitons  are  emitted 
from  strongly  excited  nonlinear  waveguides’  In  1986.  Wright  et  al.  showed  that  a 
single  waveguide  with  a  nonlinear  cover  medium  emitted  one  or  more  spatial  solitons 
into  the  cover  medium  when  excited  loo  strongly.  This  spatial  soliton  travels  away 
from  the  waveguide  at  an  angle  which  increases  with  increasing  incident  power. 
The  question  posed  then  was  whether  a  second,  parallel  waveguide  could  trap  a 
reasonable  fraction  of  the  spatial  soliton  emitted  from  the  first  guide.  This  possibility 
was  the  impetus  for  the  research  described  in  this  section.  Clearly,  as  a  minimum 
requirement,  the  spacer  region  between  two  guiding  films  must  contain  a  self- 
focusing  nonlinearity  for  spatial  solitons  to  exist  between  the  two  waveguides. 
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The  approach  used  is  the  Beam  Propagation  Method  (BPM).  Each  linear  guide 
in  the  geometry  shown  in  Fig.  7  is  defined  by  a  refractive  index  step  An 
superimposed  on  a  background  substrate  index  no-  This  involves  solving  numerically 
for  the  total  field  envelope  E(x.z)  using  the  slowly  varying  phase  and  amplitude 
equation10 


a2 

d? 


a2 

d? 


+  2ino£ 


(n*(x.|  E(x,z)| 5) 


E(x.z)  -  0 


(7) 


where  all  of  the  spatial  dimensions  are  scaled  to  k^u/c).  The  refractive  index 
distribution  n*(x.|  E(x,z)| 2)  consists  of  two  parts:  the  first  is  the  linear  index 
distribution  as  described  above,  and  the  second  is  the  nonlinear  part,  which  is 
written  as  a(x)|  E(x.z)| 2.  Equation  7  was  solved  numerically  using  up  to  40% 
transverse  grid  points.10 
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Fig.  7.  Geometry  oj  the  soliton  coupler.  Normalized  parameter  values  used  are  d  =  16, 
D  -  22,  n0  s  J.55 ,  fn  *•  0.02,  and  a *  *  1.  Guide  1  is  the  input  guide. 


Using  an  effective-particle-theory  approach,  the  conditions  on  a(x)  for  effective 
trapping  of  the  solitons  in  the  second  waveguide  were  deduced.  A  useful  geometry 
is  to  have  self-focusing  nonlinearities  both  in  the  region  between  the  two  waveguides, 
and  in  the  second  film.  Figure  8  shows  the  transmitted  flux  through  the  input  guide 
(guide  1,  Figure  8a)  and  the  parallel  guide  (guide  2.  Figure  8b)  as  a  function  of  input 
flux  into  guide  1.  Transmission  of  100%  in  guide  1  occurs,  followed  by  a  sharp 
threshold  beyond  which  soliton  emission  from  guide  1  occurs.  In  comparison,  no 
energy  is  transmitted  into  guide  2  below  the  critical  flux  for  emission  of  solitons 
from  guide  1.  but  just  above  a  significant  fraction  (a80%)  of  the  input  energy  is 
transferred  to  guide  2.  Figure  8  clearly  shows  an  extremely  sharp  switching 
characteristic  which  is  due  to  the  fact  that  soliton  emission  is  a  threshold  effect.  It 
is  also  interesting  to  note  that  the  energy  transferred  to  guide  2  stays  fairly  constant 
above  the  switch  point,  indicating  some  limiting  action. 
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Input  Flux 

Fig.  8.  Transmitted  flux  versus  input  flux  for  (a)  guide  1,  and  (b)  guide  2. 


The  phase  characteristics  of  the  transferred  guided  wave  have  not  yet  been 
analyzed.  However,  it  appears  that  the  system  may  be  quite  noisy.  For  example,  the 
BPM  results  for  the  energy  as  a  function  of  transverse  and  propagation  coordinates 
are  shown  in  Fig.  9.  Immediately  after  transfer  between  the  two  waveguides,  the 
field  distributions  undergo  oscillations  which  do  not  damp  out  for  many  hundreds  of 
wavelengths  of  propagation. 

SUMMARY 

We  have  examined  theoretically  two  nonlinear  phenomena  with  the  potential  for 
locking  together  two  waveguides/lasers. 

Clearly,  phase  conjugation  of  a  plane-wave  incident  field  onto  a  waveguide 
supporting  two  guided  pump  beams  is  a  potentially  useful  phenomenon.  This  process 
can  be  very  power  efficient  because  of  the  strong  spatial  confinement  of  the  pump 
beams.  Multi-reflection  effects  within  the  guiding  film  lead  to  oscillations  in  the 
signal  with  increasing  waveguide  thickness  and  incidence  angle. 

Coupling  between  two  waveguides  can  be  achieved  by  the  exchange  of  spatial 
solitons  when  a  self-focusing  medium  is  used  between  two  parallel  waveguides. 
Model  calculations  based  on  the  powerful  Beam  Propagation  Method  have  indicated 
transfer  efficiencies  as  large  as  80%. 
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2.25.  The 


Fig.  9.  Evolution  of  the  intensity  profile  in  a  soliton  coupler  for  5in  * 
longitudinal  coordinate  is  z  in  units  of  free- space  wavelengths. 
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